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Abstract 

We present thermal expansion a, magnetostriction and specific heat C measurements of TICUCI3, which shows a quantum phase 
transition from a spin-gap phase to a Neel-ordered ground state as a function of magnetic field around Hco — 4.8 T. Using 
Ehrenfest's relation, we find huge pressure dependencies of the spin gap for uniaxial as well as for hydrostatic pressure. For T — > 
and H ~ Hco we observe a diverging Griineisen parameter F(T) = ct/C, in qualitative agreement with theoretical predictions. 
However, the predicted individual temperature dependencies a(T) and C(T) are not reproduced by our experimental data. 
© 2008 Elsevier B.V. All rights reserved. 
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1. Introduction 

Low-dimensional quantum magnets show very rich and 
fascinating physical properties [1]. As a starting point, one 
may consider isolated spin-1/2 dimers with an antiferro- 
magnetic coupling J causing a singlet ground state and an 
excited triplet state, which are separated by an energy gap 
A = J. If such dimers are magnetically coupled to each 
other, a multitude of different theoretical models can be 
constructed, depending on the strength and the geometric 
arrangement of the interdimer coupling(s) J'. As a conse- 
quence of one (or more) non-zero J', finite dispersion(s) of 
■ the triplet state evolve along the respective direction(s) in 
reciprocal space, i.e., the excited triplets may hop along dif- 
ferent directions. Of particular interest are one-dimensional 
(ID) chains with alternating couplings J and J' between 
neighboring spins, because of qualitatively different excita- 
tion spectra of the alternating (J' ^ J) and the uniform 
(J 7 = J) chain [2]. Another example of ID coupled spin 
dimers is represented by so-called spin-ladders with the 
couplings J± and Jj| along the rungs and legs of the lad- 
der, respectively. The excitation gap of two-leg ladders is 
finite, while it vanishes for a three-leg spin-ladder [3]. This 
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difference is easily understood in the limit J± S> J\\ , where 
the two-leg ladder can be viewed as weakly coupled dimers 
and the three-leg ladder as an effective 5=1/2 chain with 
uniform chain coupling Jn . With increasing number of legs, 
the spin ladders approach the 2D antiferromagnetic square 
lattice. Another well-studied system of 2D-coupled spin 
dimers is the 2D Shastry-Sutherland model [4], which can 
be generated from the 2D square lattice by introducing one 
additional diagonal coupling Jd on every second square. 
The triangular arrangement of one Jd and two J causes a 
strong frustration and for J/ Jd < 0.7 the product state of 
singlets on every diagonal is the exact ground state. 

The above-mentioned models have been studied very in- 
tensively by theoretical as well by experimental physicists 
during the last decades, since a large number of materials 
became available which rather well represent various types 
of these models [1]. Some cuprate examples are: the 2D 
square lattice realized by the parent compound La2Cu04 
of the High-T c 's, the spin-Peierls system CuGe03 [5], spin- 
1/2 chain and ladder compounds, such as Sr2Cu03 and 
SrCu203, respectively, as well as Sri4Cu2404i containing 
both, spin chains and spin ladders [6], or the 2D Shastry- 
Sutherland model which is realized in SrCu2 (1303)2 [7]. 

TICUCI3 is another quantum spin system which has been 
intensively studied in recent years [9,10]. From the struc- 
tural point of view this compound contains a ladder-like 
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Fig. 1. Left: Phase boundary of TICUCI3 obtained from magne- 
tostriction (•) and thermal expansion (o). The dashed line is a fit 
of the form H C (T) = H C o + a ■ T<>> for T < 2 K, which yields 
Hco = 4.82 T and rf> = 1.85, while the solid line is a fit up to 8 K 
yielding <f> = 2.1. The T dependence of <f> is shown in the inset by the 
plot of Hq — Hco versus T on double-logarithmic scales, and is also 
obtained by Quantum Monte Carlo calculations [8]. Right: Schematic 
view of the phase diagram around the quantum phase transition 
from a spin-gap to a Neel-ordered ground state at Hco- The dotted 
line shows the spin-gap closing due to the Zeeman splitting, and the 
dashed lines indicate the region of enhanced quantum fluctuations 
at T > 0. The thick solid line is the phase boundary Tjf{H) and the 
thin solid lines indicate the region of enhanced thermal fluctuations. 

arrangement of Cu 2+ ions [9] . The main magnetic coupling 
J ~ 5.5 meV is present along the rungs, but there are 
various additional, rather large couplings J' present along 
different other lattice directions [11,12,13,14]. Thus, the 
magnetic system of TICUCI3 should be viewed as a set of 
thrcc-dimensionally coupled spin dimers. The inter-dimer 
couplings J' cause a strong dispersion of the triplet exci- 
tations, and as a consequence the minimum singlet-triplet 
A m ~ 0.7 meV is much smaller than J. A moderate field 
of about 5 T is already sufficient to close A m and induces a 
Neel order with staggered magnetization perpendicular to 
the applied field. If there is no magnetic anisotropy in the 
plane perpendicular to the applied field, this transition is 
in the same universality class as the Bose-Einstein conden- 
sation (BEC) and it is often termed a BEC of magnons. In 
the zero-temperature limit, it represents an example of a 
quantum phase transition [15], whose control parameter is 
the magnetic field strength (see Fig. 1). In the vicinity of a 
quantum critical point (QCP) anomalous temperature de- 
pendencies are expected for various physical properties, as 
e.g. specific heat C, susceptibility Xi thermal expansion a, 
(and resistivity p for metals) [16] . In particular, a divergence 
of the so-called Gruneisen parameter T — a/C is expected, 
when a pressure-dependent QCP is approached [17]. Ex- 
perimentally, a diverging T(T) has indeed been observed 
for different heavy-fermion compounds [18,19]. Since the 
phase transition of TICUCI3 is extremely sensitive to pres- 
sure [20,21,22,23] and the control parameter may be easily 
tuned by a variation of the field, this compound is ideally 
suited to study such generic properties of a QCP. 

We present high-resolution measurements of the uniaxial 
thermal expansion 014 = d In Li /dT and the magnetostric- 
tion ei = \Li(H) — Lj(0)]/-L,(0) along different lattice di- 
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rections i {Li is the respective sample length along i) as 
well as specific heat C and magnetization M data. The 
length changes have been measured down to 250 mK by 
a home-built capacitance dilatometer and C by a home- 
built calorimeter for T > 500 mK, while the magnetization 
has been studied by a commercial vibrating sample mag- 
netometer (Quantum Design) for T > 1.9 K. All proper- 
ties have been studied in magnetic fields up to 14 T. Since 
TICUCI3 easily cleaves along the (010) and (102) planes of 
the monoclinic structure, we measured Li{H,T) perpen- 
dicular to these planes on a single crystal of dimensions 
1.7 x 1.5 mm 2 perpendicular to (010) and (102), respec- 
tively. In addition, the [201] direction, which is perpendic- 
ular to both other directions, was measured on a second 
crystal of length Lpoi] = 4.4 mm. For all three measure- 
ment directions i the magnetic field was applied along the 
same direction, namely perpendicular to the (102) plane. 

2. Results and Discussion 

Fig. 2 shows cti measured along all three directions for dif- 
ferent magnetic fields. In zero field, cti has n0 visible anoma- 
lies, but is strongly anisotropic. For H > 5 T, pronounced 
anomalies develop and shift systematically towards higher 
T with increasing H. The ccj curves for i = (010) and (102) 
agree well with our previous results measured on a different 
crystal for T > 3 K [22,24]. The anomalies of a, signal large 
uniaxial pressure dependencies of TV, which are related to 
a and C via Ehrenfest's relations 

dT N Aa, dH c A\ 

— — = V m T N —— and — — = V m — — . (1) 

opt AC dpi A\ 

Here, V m is the molar volume and Ao^ and AC denote the 
respective mean- field jumps at Tn ■ The second expression 
of Eq. (1) relates the pressure dependencies of the transi- 
tion field to the jumps of Xi = dei/dH and of the differen- 
tial magnetic susceptibility x — dM /dH. For i — (010) and 
(102), the dTjy/dpi largely cancel each other under hydro- 
static pressure, since the anomalies are of comparable mag- 
nitudes but of opposite signs. Thus, the hydrostatic pres- 
sure dependence dTjqjdphyd.ro is essentially determined by 
the sign and size of the anomaly of 0:201- Obviously, the 
anomalies of 0:201 are the largest ones and their positive 
signs mean that T/v drastically increases for uniaxial pres- 
sure along [201] as well as for hydrostatic pressure. Fig. 3 
shows the field derivatives of the magnetostriction. Again, 
we find the characteristic anisotropy that the anomalies of 
A010 and A 10 2 are of similar sizes but opposite signs, while 
significantly larger anomalies are present for A2oi- Thus, 
the hydrostatic pressure dependence of He is again essen- 
tially identical to that for uniaxial pressure along [201]. 

In Ref . [22] we have shown that the pressure-dependencies 
of T/v, He and the magnetic susceptibility in the paramag- 
netic phase can be traced back to the pressure-dependent 
changes of a single parameter, which in the case of TICUCI3 
is the intra-dimer coupling J. This conclusion was based 
on the measurements perpendicular to the (010) and (102) 
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Fig. 2. Thermal expansion along different directions for various mag- 
netic fields applied perpendicular to the (102) cleavage plane of 
TlCuCl 3 . Note the different scale for the [201] direction. 

planes, and is fully confirmed by the additional new data 
along the [20 f] direction. From the zero-temperature ex- 
trapolations of the uniaxial pressure dependencies of He 
for i = (010) and (102) we estimated <91n A m /dpi ~ 
+ 190 %/GPa and ~ -180 %/GPa. Since the anomalies 
for the [201] direction are about twice as large, we obtain 
9 In A m /dpi ~ —370 %/GPa for pressure along [201] and 
~ —360 %/GPa for hydrostatic pressure. The latter value 
is in reasonable agreement with direct measurements un- 
der hydrostatic pressure, which yield 9 In A m /dph v dro — 
—400 %/GPa for the initial slope at ambient pressure [23]. 

The shape of the aj anomalies is typical for a second- 
order phase transition with a pronounced mean-field con- 
tribution causing a jump Aa at TV, superimposed by fluc- 
tuations causing a divergence a oc t v with the reduced 
temperature t = \T — T^\/T^ and the critical exponents 
v depending on the universality class of the phase transi- 
tion. On approaching H C q = H C (T -> 0) ~ 4.8 T, the 
cx,i anomalies broaden to some extent (see below). The Aj 
anomalies also show a pronounced fluctuation contribution 
for T > 2 K, but become more jump-like for lower T. The 
changing shapes of both, the Aj and the a, anomalies can be 
intuitively understood from Fig. 1, because (i) the absolute 
temperature region around the phase boundary where fluc- 
tuations become important decreases with decreasing TV, 
and (ii) close to Hco the phase boundary is crossed with a 
very small slope as a function of T, so that the fact that He 
is not infinitely sharp becomes more and more important. 
As a criterion for He, we used the maximum of the sec- 
ond derivatives d 2 €i/dH 2 , whose full widths at half maxi- 
mum amount to ~ 0.45 T (see Inset of Fig. 3). We suspect 
that this width mainly arises from internal stresses, which 
broaden the transition due to the strong (uniaxial) pressure 
dependencies of He- We have also investigated, whether 
there is a finite hysteresis of He by measuring (H) with 
increasing and decreasing H. For a drift rate of ±0.1 T/min 
we obtain a difference H c p - H c own ~ 0.04 T, which does 
hardly change with temperature. At T = 0.3 K we also 
measured with ±0.01 T/min and found H c p - H c own ~ 
0.01 T, i.e. , the observed hysteresis partly arises from the 
finite field drift. Thus, we regard the phase transition as 
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Fig. 3. Magnetostriction coefficient Xidti/dH measured along differ- 
ent directions i at various temperatures. In all cases the magnetic 
field has been applied perpendicular to the (102) plane (see text). 

a second-order one with a weak first-order contribution. 
The first-order contribution most probably arises from the 
large spin-lattice coupling, which may drive a second-order 
into a first-order transition [25]. It is also possible that a 
transition transforms from second to first order, when T c 
is suppressed towards K by an external parameter. How- 
ever, the weak temperature-dependence of the hysteresis 
observed in TICUCI3 does not give any evidence for such 
a scenario, and we suspect that the transition of TICUCI3 
remains (almost) continuous down to lowest T. 

As mentioned above, a diverging Griineisen parameter 
T(T -> 0) has been predicted at H = H C a [17]. Before 
comparing our experimental data to this prediction, we will 
discuss from a phenomenological point of view what may 
be expected for T when the QCP is approached along dif- 
ferent routes in the phase diagram. Using Maxwell's rela- 
tions one finds that a = -dS/dp, while C/T = dS/dT. If a 
thermodynamic system can be described by a single energy 
scale E, its entropy S only depends on the ratio T/E, i.e. 
S(T, E) = S{T/E) with a model-dependent function S{x). 
Comparing the T- and p-derivatives of S(T/E(p)) yields 
the Griineisen scaling 



r = 



a(T) _dlnE 
C(T) ~ dp 



(2) 



which is temperature independent. Prominent examples of 
such systems are the Debye model with E = Qd, the almost 
free electron gas with E = Ep, or magnetically ordered 
states with exchange coupling J = E (for T <C T c ). In real 
systems, one usually observes a weakly T-dependent T(T), 
which is due to the fact that the above-mentioned single- 
parameter models only consider the leading energy scale 
and neglect others. If several energy scales Ei are equally 
important, the individual Ei usually have different pressure 
dependencies, and therefore the p- and T-derivatives are, 
in general, not (almost) proportional to each other. How- 
ever, it is possible that for different temperature regions 
different T^'s are dominant and in the respective regions 
r ~ dhiEi/dp holds. An example is a coupled spin-dimer 
system as it is realized in TICUCI3. For high temperatures 
(T 3> J, J'), the behavior is determined by the average 
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Fig. 4. Left: Difference of the specific heat C in finite and zero 
magnetic field. The inset shows the anomaly heights as a function of 
field. Right: Expanded view of the low-temperature C/T for different 
fields below and above the critical field Hco — 4.8 T. The arrow 
indicates increasing magnetic-field strength, (see text) 
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Fig. 5. Magnetization (symbols) of TICUCI3 at different tempera- 
tures. The finite magnetization in the gap region can be well repro- 
duced by Brillouin functions (lines) assuming magnetic impurities 
with spin 5" = f . (see text) 



spin gap (A), while the minimum gap A m becomes dom- 
inant at T <C A m . Thus, T(T) varies from ~ d\nA m /dp 
to ~ d\n(A)/dp with increasing T. In TICUCI3, the zero- 
field gap A°„ ~ 8 K linearly decreases with H [13], i.e. 
A m (H) = A^ — h where h — g H/Jcb, g is the ^-factor and 
k,B Boltzmann's constant. Thus, at a given field H < Hco 
the temperature, below which T is expected to approach 
a constant, decreases with H. For H — ► Hco, however, 
T(T < A m ) ~ A J"( g) -Q^r diverges, because A m -> 0. 

For H > Hco, i-e- in the ordered phase, the characteris- 
tic low-temperature energy scale is given by the spin wave 
velocity v. To our knowledge, the exact dependence of v 
on H — Hco is not known, but it is quite natural that v 
disappears when Hco is approached. For simplicity, we as- 
sume v ex (H — Hco) n with n > and derive T(T — > 0) ~ 
d In v/dp oc H _ 1 Hcg ~ d g p Q0 ■ Thus, the Griineisen parameter 
is again expected to approach a constant, which diverges 
for H — > Hco , but the sign of the divergence for H > Hco 
is opposite to that for H < Hco- 

In a next step we will approach the QCP along the phase 
boundary Tjv(-ff) for H > Hco- For clarity, we use the 
approximation TV (if) — b(H — Hco) v with ip = 1/0 — > 
2/3 for H — > Hco (see Fig. 1 and Refs. [8,26]), and calcu- 
late dT N /dp = -ipb (H - Hco) {ip ~ 1] ■dH C o/dp. Obviously, 
dTjsi/dp diverges for H — > Hco for <p < 1, and the sign of 
this divergence is opposite to the sign of dHco / dp. We em- 
phasize that this result follows from the infinite slope of the 
phase boundary for H — > Hco, an d does not depend on the 
particular choice of T^(H). Since the pressure dependence 
oiTpj is given by Ehrenfest's relation (1), the ratio Aa/AC 
of the thermal-expansion and specific- heat anomalies at TV 
is expected to diverge for H — > Hco- From Eq. (1), it is 
also clear that the vanishing TV would cause a divergence 
of Aa/AC even if the slope of the phase boundary was fi- 
nite and <9T/v /dp would thus not diverge for H — > Hco- 

Let us summarize the above considerations. On ap- 
proaching Hco w e expect (i) for H < Hco that T(T <C 
A(7T)) -> dhxA(H)/dp, which diverges for IT -> H C o, (n) 
for IT > Hco a similar divergence of T(T — > 0), but of the 



opposite sign, and (iii) a divergence of the ratio Aa/AC, 
which has the same (opposite) sign as the divergence of 
r above (below) Hco- Since the above considerations are 
rather general, one may expect a divergence of T(T) close 
to many kinds of transitions, whose T c is suppressed to 
K. To obtain more information about a quantum phase 
transition, one has to consider the actual temperature de- 
pendencies a(T), C(T), and/or T(T), as it has been done 
e.g. by the authors of Refs. [17,27,28]. 

In Fig. 4 we show the specific heat anomalies for differ- 
ent magnetic fields. In agreement with Ref. [29] we find 
rather small anomalies even for the largest field, and their 
magnitude rapidly decreases when Hco is approached (see 
Inset of Fig. 4). Since the magnitude of the respective a« 
anomalies changes much less with field (see Fig. 2), the ex- 
pected divergence of dT^/dpi for H — > Hco is obviously 
confirmed by the experimental data, since the denomina- 
tor in Eq. (1) vanishes. This is the case for all three di- 
rections of uniaxial as well as for hydrostatic pressure. In 
the right panel of Fig. 4 we show an expanded view of the 
low-temperature behavior of some C(T, H)/T curves. For 
zero field, the onset of an anomaly can be clearly seen [30], 
and this anomaly is suppressed above about 3 T. We sus- 
pect that this anomaly arises from an ordering of magnetic 
impurities. The presence of such impurities is also evident 
from the finite magnetization in the gap region H < Hco 
at low T, which can be well reproduced by Brillouin func- 
tions (see Fig. 5). The corresponding fits of the data at T = 
1.9 K and 2.5 K yield 0.4 % of magnetic impurities with spin 
S = 1. The Brillouin function calculated for the same pa- 
rameters and T = 4 K is somewhat smaller than the exper- 
imental data for H > 1 T. This is expected, because at this 
higher temperature a sizeable magnetization from excited 
triplets is already present. Probably, the S = 1 impurities 
are mostly ferromagnetically aligned spin dimers, because 
the intra-dimer coupling between the spin-1/2 Cu 2+ ions 
arises from a ~ 96° superexchange via the p orbitals of the 
Cl _ ions, which is very sensitive to changes in the bond an- 
gle. The impurities strongly influence the low-temperature 
behavior of C(T) for low fields, but become much less in- 
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Fig. 6. Top: Thermal expansion perpendicular to the (010) plane for 
different magnetic fields below (left) and above (right) the critical 
field. The arrows indicate increasing magnetic-field strength. The 
bottom panels show some cti/T curves for the other directions. 

fluential at higher fields because the moments are almost 
completely saturated for H > 3 T and T < 2 K. It is, how- 
ever, unclear to what extent the impurities may change the 
(critical) behavior close to the phase transition. 

Fig. 6 shows an expanded view of the low-T behavior of 
oti/T. In zero field, cti/T continuously approaches zero for 
T —> 0, while for larger H it shows a pronounced shoul- 
der, which systematically increases and reaches a maximum 
slightly below Hco- For larger fields clear anomalies with a 
sign change of oti/T occur, and these anomalies systemati- 
cally sharpen with further increasing field. The behavior is 
essentially the same for all three directions, only the magni- 
tudes and signs are different [31]. As already mentioned, we 
attribute the broadening of the anomalies when Hco is ap- 
proached from larger fields to the finite width of the phase 
transition. This also explains that the ai(T) curves show 
anomalies already for H > 4.6 T, i.e. below Hco — 4.8 T 
determined by the magnetostriction measurements. 

In Fig. 7 we present Ti for different magnetic fields. For 
all three directions Ti shows the tendency to diverge with 
decreasing T for H ~ 4.5 T. For lower H, the r,(T) follow 
the same curve at higher T, but seem to approach finite 
values for T — > 0. The magnitudes of these limiting values 
increases with increasing field and the temperature, below 
which the deviation sets in, decreases. The Ti(T) for H > 
Hco also follow the general curve at higher T, until a large 
anomaly signals the crossing of the phase boundary. For all 
three directions, the magnitudes of these anomalies drasti- 
cally increase with decreasing H and the signs are opposite 
to the respective signs of the diverging Ti(T) for H < Hco- 

On this qualitative level, our experimental data of Ti(T) 
very well confirm the behavior, which one can expect from 
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Fig. 7. Griineisen parameters Fi obtained from the thermal expansion 
along different measurement directions i in magnetic fields below 
(symbols, increasing field strength is indicated by the arrows) and 
above (lines) the QCP at Hco — 4-8 T. The upper right panel shows 
Ti(T) for H < Hco on double-logarithmic scales. Here, the negative 
Ti(T) for t = (102) and [201] have been divided by a factor of -2 
and —3, respectively. The solid line is a power-law fit of Toio yielding 
I/T 4 / 3 and, for comparison, the expected 1/T behavior is shown by 
the dashed line (see text). 

the above considerations of the gap closing for H < Hco 
and the shape of the phase boundary for H > Hco- For a 
deeper understanding, one has to compare the experimen- 
tal data quantitatively to theoretical predictions. Accord- 
ing to Ref. [17] the following temperature dependencies are 
expected H = Hco'- 



C/TocVT, a/T oc 1/Vf , and Tex 1/T 



(3) 



In the upper right panel of Fig. 7 we show the diverging 
Ti(T) for H < Hco on double- logarithmic scales. Because 
of the negative signs, r i0 2 and T2oi have been divided by 
the factors —2 and —3, respectively. Apparently, the slope 
is the same within experimental accuracy. The solid line is 
a power-law fit of Toio which yields T -4 / 3 and describes the 
experimental data reasonably well for about one decade. 
For comparison, the predicted 1 /T behavior is also shown 
(dashed line). In view of the fact that the theoretical pre- 
diction only considers the irregular contributions of C and 
ctj, while the experimental data also contain the phononic 
contributions of C and aj , one may tend to the conclusion 
that our data nicely confirm the theoretical prediction for 
r(T). However, the agreement between theory and exper- 
iment becomes much worse when the individual tempera- 
ture dependencies of C/T and Oj/T are considered. Nei- 
ther the data of Fig. 4 nor those of Fig. 6 give any indi- 
cation to follow the predicted temperature dependencies 
of Eq. (3). Concerning the specific heat data, one might 
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argue that the predicted VT behavior is difficult to iden- 
tify because of the phononic contribution and the influence 
of the magnetic impurities. This argument is less convinc- 
ing for a.i(T)/T, since (i) the predicted divergence should 
be seen despite a (regular) phononic contribution and (ii) 
the ordering of the magnetic impurities does not cause a 
sizeable anomaly in the zero-field data [31]. Thus, we con- 
clude that our present experimental data of TICUCI3 do 
not confirm the theoretical predictions [17]. However, our 
data do not disprove the theoretical predictions either. Ex- 
perimentally, one can suspect that in order to observe the 
predicted temperature dependencies it would be necessary 
to study (i) samples with significantly reduced transition 
widths and less magnetic impurities, and (ii) it might be 
necessary to extend the measurements to lower tempera- 
tures. Moreover, the theoretical predictions have been cal- 
culated for clean and isotropic systems. It is not clear to 
what extent the temperature dependencies of a and C are 
influenced by disorder or a finite magnetic anisotropy. The 
latter is reflected in the ~ 10 % anisotropy of the g factors 
of TICUCI3 for different magnetic field orientations [32]. 

3. Summary 

In summary, we have presented high-resolution measure- 
ments of thermal expansion and magnetostriction along dif- 
ferent lattice directions of TICUCI3 . Both quantities show 
very pronounced and strongly anisotropic anomalies at the 
phase boundary of the field-induced Neel order Tn(H) for 
H > Hco, and signal very large and strongly anisotropic 
uniaxial pressure dependencies of the transition tempera- 
tures and fields. The hydrostatic pressure dependence of 
the spin gap A m for H — > Hco obtained from our data us- 
ing Ehrenfest's relations is in reasonable agreement with 
the value observed by direct measurements under hydro- 
static pressure. In addition, our data confirm the diverg- 
ing pressure dependencies dTw/dpi for H — > Hco, which 
are expected from the infinite slope 8Tn /dH of the phase 
boundary for H — ► Hco, i- e - when the quantum critical 
point is approached. For H < Hco, the Gruneisen param- 
eters Ti — oti/C are expected to approach constant values 
for T < A(H), which diverge for H -> H co , and Ti(T) oc 
1/T has been predicted at H = Hco- In fact, the exper- 
imental Ti(T) for all three measurement directions are in 
qualitative agreement with these expectations. However, 
the temperature dependencies predicted for the individual 
quantities cti(T) and C(T) are not observed experimentally. 
For H ~ Hco the low-T behavior of both a, and C is in- 
fluenced by the finite transition width and for lower fields 
at least C(T) is also affected by the presence of magnetic 
impurities. Thus, future measurements on samples of im- 
proved quality as well as calculations considering the influ- 
ence of disorder and weak magnetic anisotropy may clarify 
the reasons for the puzzling temperature dependencies of 
oh(T), C(T), and T^T). 
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